Fluctuation theorem for black-body radiation 
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The fluctuation theorem is verified for black-body radiation, provided the bunching of photons is 
taken into account appropriately. 
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I. INTRODUCTION 

The study of black-body radiation has played a promi- 
nent role in the discovery of quantum mechanics. More 
recently, the quantum features of light have been revealed 
in many beautifuU optical experiments, and found to be 
in full agreement with the quantum theory of electro- 
magnetic radiation [TJ [SJ [3]. In particular, we cite the 
bunching of photons. Bunching refers to the tendency of 
photons to arrive together at closely spaced detectors, re- 
sulting in super-Poissonian counting statistics (variance 
large than the mean). While a classical explanation is 
available Hj, the detailed understanding and description 
of the phenomenon requires quantum field theory [U [5] . 
The bunching of photons, and more generally of bosons 
[B], is usually explained in terms of the more familiar 
quantum statistical exchange force which introduces an 
effective attraction between identical bosons, when their 
wavefunctions overlap [7] . The phenomenon of bunching 
was first observed by Hanbury Brown and Twiss 8 as 
the correlation between intensities of light from a single 
star falling on different detectors. Bunching of photons 
in a single detector experiment was later observed [9 
and found to be in agreement with a simplified quantum 
mechanical argument by Glauber [TU]. For the case of 
equilibrium black-body radiation, the photon count has 
been calculated explicitly [5] and is given by a negative 
binomial distribution rather than the Poissonian distri- 
bution which describes independent arrivals. 

In this Letter, we show that this result for the photon 
counting statistics is essential to find agreement with a 
recent result from nonequilibrium statistical mechanics, 
namely, the so-called fluctuation theorem. This theo- 
rem states that the probability distribution P{/S.S) to 
observe an entropy production AS* during a time inter- 
val t in a nonequilibrium steady state obeys the follow- 
ing symmetry relation for asymptotically long values of 

t [n nu [n El [13 [n ini : 



_P(ASO_ 
P(-AS') 



' exp{AS'/fc}. 



(1) 



In words, the probability of observing a positive entropy 
change is exponentially larger than that of the corre- 
sponding negative change. When the system starts in 
a state of canonical equilibrium which is subsequently 
perturbed by a time-dependent change of the Hamilto- 
nian, the above fluctuation theorem, referred to as the 
transient fluctuation theorem, is valid for all times, and 



not just asymptotically large ones [T^ITH] . 
The fluctuation theorem has been verifled in several the- 
oretical [Il[ini[lTllll[23[2lll5l[2SJ[13!^/^/3S[n^ 

and experimental settings [SlllMllMl [35J E6 , 37, 38J. The 
physical origin of the fluctuation theorem is to be found 
in the time-reversal symmetry of the underlying Hamil- 
tonian dynamics. As such, this result can be viewed as a 
generalization of the ideas of Onsager [55 ]. 
We shall investigate the steady state of radiation ex- 
change between two black bodies, at equilibrium at differ- 
ent temperatures, when they are connected to each other 
through a small aperture. To present the problem and 
to contrast the role of the quantum mechanical bunching 
for photons with that of classical particle exchange, we 
start by studying the classical counterpart of effusion of 
the ideal gases. 



II. EFFUSION OF A CLASSICAL IDEAL GAS 

Consider two infinitely large reservoirs A and B, each 
containing classical ideal gases at (canonical) equilibrium 
with temperatures Ta and Tb and densities pA and ps, 
respectively. We perform the following experiment (see 
fig. [l]). During a fixed time interval t, a small hole of 
surface area a between the reservoirs is opened. We as- 
sume, for better comparison with the photon crossings, 
that the opening contains an energy filter, allowing only 
particles with kinetic energy in the range Eq ± 5E/2 to 
move across the hole. The net amount of energy /S.U and 
net number of particles A7V transferred from AtoB are 
measured. We consider the limit of a small energy win- 
dow, SE Eq, so that AU = EqAN, and the exchange 
in the number of particles is the only relevant variable. 
Since the hole is small enough so that the canonical equi- 
libria in the reservoirs are not significantly perturbed, the 
corresponding entropy change AS* is given by standard 
thermodynamics : 



AS-: 
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T^ 



1 

T^ 



Eq 



g-ghA^. (2) 



Inserting the well-known espressions |40j for the chemical 
potentials fia of the ideal gases in reservoirs a E {A, B}, 
one finds: 
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FIG. 1: (Color online) Set-up for the effusion of ideal gases: 
a) Two separate ideal gases in equilibrium at their respec- 
tive temperature and density, b) During a time t, a small 
hole is opened and the net transfer of energy AC/ and AA'' is 
measured. 



Since the particle exchange AA^ is the result of individ- 
ual gas particle crossings, it is obviously a random vari- 
able, and hence so is AS*. The probability distribution 
P(AiV) can be calculated as follows. We define Pa{j) 
as the probability that j particles leave the reservoir a 
in the specified time interval. As all the particles move 
independently of one another, Pa{j) is Poissonian |41j : 



(4) 



The escape rate can be calculated following standard 
arguments from kinetic theory (see fig. [2]). We get 



v^^^i^^e-^'^'^^-SE. (5) 



Since the escape of particles from the left and right reser- 
voirs are independent phenomena, P(AN) is given by the 
sum 



PiAN) = Y.pa{AN + i)pb{i)- 



(6) 



i=0 



Combining eqs. ^ and ^ leads to the following explicit 
expression for P{AN): 



P{AN) = e 



Ian {2t^vAVB) ■ (7) 



Recalling that the modified Bessel function /atv is an 
even function of AA'^, the fiuctuation theorem is verified 
by direct substitution of this result in eq. (jlj. 
We mention a simpler proof, more clearly related to the 
underlying micro-reversibility. Consider the exchange 
process, with AN + i particles going from A ^ B, and 
i particles from B ^ A. This occurs with probability 
Pa{AN + i)pB{i)- The time reversed process, with i par- 
ticles going from A ^ B and AA^ -I- i particles from 
B ^ A, has probability pA{i)pB{AN + i). Their ratio 




Vg cos(9) dt 



FIG. 2: The escape rate is determined by counting the parti- 
cles that have the correct kinetic energy and which are able 
to reach the opening in a short time dt (these are located in 
the cylinder with volume avo cos{9)dt, vq = ^ 2Eo / m) . In- 
serting the Maxwellian velocity distribution and adding up all 
contributions from the different angles leads to the expression 
eq. (|5|. 



satisfies the following detailed fluctuation theorem: 
Pa{AN -f i)pB{i) _ ( VA^ 
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using eqs. ([s]) and (|2|. The fluctuation theorem itself 
follows immediately. We have 

00 

P{AN) = ^p^(A7V + i)pBW 



i=0 

00 / 

i=0 ^ 



AN 



PA{i)pB{AN + i) 



VA 
VB 



AN 



P(-AA). 



(9) 



III. BLACK-BODY RADIATION 

We now turn to a similar set-up for the study of fluctu- 
ations in black-body radiation. Two large empty cavities 
{A and B), whose walls are kept at a fixed, but different 
temperatures, Ta and Tg respectively, act as sources of 
black-body radiation (see fig. |3|. During a fixed time in- 
terval t, radiative exchange becomes possible by opening 
a small aperture (surface area a) , which permits free pas- 
sage of photons with a frequency in the range ujq ± (50^/2. 
We again consider the limit of a monochromatic filter 
5lj <^ luq. Following standard thermodynamics, the en- 
tropy change AS*, upon transfer of a net number of pho- 
tons AA^ from A to B (and hence of energy fiwoAN), is 
given by: 



AS = 



1 



1 



hwoAN. 



(10) 
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FIG. 3: (Color online) Set-up for black-body radiation: a) 
Two cavities are equilibrium each at its own temperature, b) 
During a time t, a small hole is opened and the net transfer 
of energy AU is measured. 



Photons coming from the two different reservoirs are in- 
dependent of each other. Hence the probability P{AN) 
to observe a net transfer of AA^ photons from A to i? is 
again given by eq. ([6|. However, photons coming from 
the same reservoir are not independent, and their escape 
is no longer governed by a Poisson distribution. The cal- 
culation of Pa{j), which is also referred to as the photon 
counting distribution, requires a fully quantum mechan- 
ical description of the electromagnetic field [Ij |5]. For 
large times {tSu ^ 1), it is given by the following nega- 
tive binomial distribution: 

P.ij) = (1 - e-^-"'^-^y' (e-^-o^^-^y , (11) 

where v is defined as 



(27rc)2 ■ 



(12) 



It is now a matter of simple algebra to verify that the 
following detailed fluctuation theorem is obeyed: 



PA(AiV-Hi)pj3(i) 
PA(«)ps(AA^ + i) 



J--^)ri(.c.oA7V/fc 



(13) 



Consequently, following the steps of eq. ([9]) , the fluctua- 
tion theorem is also verified. 

We emphasize that the correlations between the pho- 
tons, resulting in the negative binomial distribution, is 
an essential ingredient. A semiclassical approach, which 
assumes statistical independence between photons, can- 
not be reconciled with the fluctuation theorem (see Ap- 
pendix) . 



IV. DISCUSSION 

We close with a number of comments on the above re- 
sults. 

We have presented the analysis above for the case of a 
monochromatic window connecting the reservoirs. This 
permits the derivation of a detailed fluctuation theorem. 
It is also revealing to note that the entropy production 
becomes identically zero for the classical gas for a par- 
ticular choice of the energy window: the r.h.s. of eq. ^% 



is identically equal to unity for a particular value of the 
energy, namely. 
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At this specific value of the kinetic energy, the one- 
particle energy distributions in reservoirs A and B cross 
each other, so that the particles with this energy have 
the same number density on both sides. Such an equilib- 
rium state can be used to reach optimal thermodynamic 
Carnot and Curzon-Ahlborn efficiencies [H]. This type 
of equilibrium can also be achieved with electrons [43], 
but not with photons because the latter have zero chem- 
ical potential. 

The fluctuation theorem for the effusion of an ideal gas 
is valid for all times, while for photons the theorem has 
been proven here only for large times. This is consistent 
with the observation that the effusion of an ideal gas is a 
process without memory. In such a case, the distinction 
between the steady state and transient versions of the 
fluctuation theorem disappears. A transient theorem for 
photons valid for all times can also be obtained in prin- 
ciple, but this would require the exact evaluation of the 
transient photon count, starting from zero and progress- 
ing to the steady state, after the aperture is opened. 
As expected from the general argument connecting 
bunching to the statistical exchange interaction, one finds 
that fermions display anti-bunching [33] . While this phe- 
nomenon has no classical analogue, it is surprising to 
find, for example, that the statistics of an electron cur- 
rent through a small channel can be fully reproduced by 
a simple classical random walk model, namely, the sym- 
metric exclusion process [31]. For this simple model, the 
fiuctuation theorem is again verified. The validity of the 
fiuctuation theorem for fermions has also been confirmed 
by quantum field theoretic calculations [45] . 



V. APPENDIX 

A semiclassical description assumes the photons to be 
independent of each other. The photon counting distri- 
bution is then given by eqs. ^ and In this case, 
the quantum aspects of the photons are only taken into 
account via their occupation number density as implied 
by Planck's law of radiation. This leads to the following 
result for the escape rate: 



{2-KcY {t 



1 



(15) 



The ratio of these rates (cf. eq. (|8|) does not yield the 
desired result expjA^/fc} with A5 given by eq. ( 10 ), and 
hence the fluctuation theorem is not satisfied. 
The picture of independent photons has been used to cal- 
culate the average number of photons leaving compart- 
ment a in a time interval t (see, e.g., [46] ) . Since this 



4 



quantity is not influenced by bunching, the correct result Integrating over 5uj, one recovers the Stefan-Boltzmann 
is obtained in this case: law. However, the higher order moments are not cor- 

rectly reproduced in the independent photon picture. 

(i) = Vat. (16) 
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